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Diakoptic  Theory  for  Multielement  Antennas 
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Abstract — A  theory  is  presented  for  the  analysis  of  multielement 
antennas  which  consist  of  interconnected,  conductive  structure 
elements  of  electrically  snull  dimensions.  The  theory  is  based  on  the 
retarded  electromagnetic  potentials  which  permit  a  diakoptic  ap¬ 
proach  to  the  problem.  The  antenna  is  broken  up  into  its  individual 
structure  elements.  Each  element  is  assumed  to  be  excited  by  currents 
which  are  impressed  at  its  terminals,  i.e.,  junctions  with  adjacent 
elements  (current  coupling)  and  by  the  electric  fields  of  the  currents 
and  charges  on  all  the  other  elements  (field  coupling).  Both  excitations 
are  treated  independently.  Each  impressed  current  produces  a 
“dominant”  current  distribution,  a  characteristic  of  the  element, 
which  can  be  readily  computed.  Current  coupling  is  formulated  by 
“intrinsic”  impedance  matrices  which  relate  the  scaler  potentials  at 
the  terminals  of  an  element,  caused  by  its  dominant  current  dis¬ 
tributions,  to  the  impressed  currents  of  the  dement.  Field  coupling 
produces  “scatter”  currents  on  all  the  elements  and  is  formulated  by 
a  “field -coupling"  matrix  which  relates  the  scalar  potentials  at  the 
terminals,  earned  by  field  coupling,  to  the  impressed  currents  at  all 
the  terminals.  Intrinsic  and  “field-coupling”  matrices  are  combined 
to  form  the  “complete"  impedance  matrix  of  the  diak opted  antenna. 
Enforcing  continuity  of  the  currents  and  equality  of  the  scalar 
potentials  at  all  the  Interconnections  between  the  elements  yields  a 
system  of  linear  equations  for  the  junction  currents  and  the  Input 
impedance  of  the  antenna.  Current  coupling  dominates  field  coupling. 
Field  coupling  is  primarily  affected  by  the  dominant  current  dis¬ 
tributions  of  the  elements,  and  in  general  the  scatter  currents  have 
negligible  effect  on  M.  Although  detailed  numerical  investigations  will 
be  presented  in  another  paper,  a  simple  example  is  included  here  to 
demonstrate  that  the  diakoptic  theory  yields  very  good  results  even  if 
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I.  INTRODUCTION 

AN  ANTENNA  which  consists  of  a  number  of  intercon¬ 
nected  conductive  structure  elements  of  electrically  small 
dimensions  is  shown  in  Fig.  1.  This  multielement  design  was 
chosen  to  obtain  a  broad-band  highly  efficient  antenna  of 
relatively  low  profile  [11. 

An  analytical  treatment  of  such  a  composite  structure  ap¬ 
pears  to  be  a  rather  hopeless  undertaking  Commonly  used 
numerical  techniques  would  require  computers  with  large 
storage  capacity  [2],  This  paper  offers  a  new  approach  to 
the  problems  of  this  kind,  which  holds  promise  for  improved 
numerical  efficiency.  According  to  this  approach  the  com¬ 
posite  structure  is  diakopted  into  its  individual  structure 
elements.  As  a  simple  example  Fig.  2  shows  a  diakopted 
dipole  with  end  capacitor  plates.  Each  structure  element  is 
characterized  by  electrical  quantities  which  depend  only  on 
the  size  and  shape  of  the  element,  and  the  assembly  is  treated 
similarly  to  the  interconnection  of  n-port  networks. 

The  excitation  of  each  element  is  ascribed  to  two  causes: 
a)  the  currents  entering  the  element  at  its  “terminals,”  i.e., 
junctions  with  adjacent  elements  or  the  source,  and  b)  the 
fields  of  the  currents  and  charges  on  all  the  other  elements. 
The  first  is  referred  to  as  “current  coupling”  and  the  second 
as  “field  coupling.”  Each  excitation  is  treated  separately. 
Current  coupling  implies  hypothetical  sources  with  a  single 
terminal  and  the  capability  of  impressing  a  current  onto  a 
conductor.  Although  such  sources  violate  the  continuity 
condition,  their  assumption  is  permissible  if  the  electromag¬ 
netic  fields  are  expressed  by  the  retarded  electromagnetic 
potentials.  Although  the  continuity  condition  is  violated  in 
the  treatment  of  individual  structure  elements,  it  is  restored 
when  the  elements  are  interconnected. 

Let  us  for  the  moment  disregard  field  coupling.  If  a  current 
is  impressed  at.  a  terminal  of  a  structure  element  the  current 
spreads  over  the  surface  of  the  element  and  produces  a  current 
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Fig.  1.  Broad-band  multielement  antenna. 


Fig.  2.  Diakopted  capacitively  loaded  dipole. 


Fig.  3.  Excitation  of  single  terminal  structure  element. 

with  a  current  1  impressed  at  the  terminal,  i.e.,  contact  area 
in  the  center  of  the  plate  (Fig.  3).  The  contact  area  a  is  con¬ 
sidered  very  small  compared  with  the  surface  of  the  element. 
Excitation  by  an  impressed  current  cannot  be  treated  with 
Maxwell’s  equations,  because  Maxwell’s  equations  imply 
sources  which  separate  positive  and  negative  charges.  In 
contrast,  impressed  currents  require  sources  which  produce 
charges.  The  regarded  electromagnetic  potentials  do  not 
impose  any  conditions  on  the  source  and  therefore,  can  be 
used  for  our  problem. 

If  i(r)  is  the  surface  current  density  and  q(f)  the  surface 
charge  density  due  to  the  impressed  current  /,  the  retarded 
potentials  are 


A(r)  = 

P 

4rr  J 

f  Kf'yG(f,r')dS 
s 

(vector  potential) 

(1) 

*'*)  = 

1 

4ire  ■ 

f  q(r')G(f,  f')dS 
's 

(scalar  potential) 

(2) 

with 


distribution  which  is  uniquely  determined  by  the  ' .  jmetry 
of  the  element  and  the  location  of  the  terminal.  There  are  as 
many  dominant  current  distributions  as  there  are  terminals. 
The  relationship  between  the  scalar  potentials  at  the  terminals, 
produced  by  the  dominant  current  distributions,  and  the 
impressed  currents  is  formulated  by  the  “intrinsic  impedance 
matrix”  of  the  element. 

Field  coupling  excites  scatter  currents  which  are  super¬ 
imposed  on  the  dominant  current  distributions.  The  scalar 
potentials  at  the  terminals  due  to  field  coupling  depend  on  all 
the  impressed  currents.  Their  relationships  with  the  impressed 
currents  are  formulated  by  a  “field-coupling”  matrix.  Intrinsic 
impedance  and  “field-coupling”  matrices,  when  combined, 
form  the  “complete  impedance  matrix"  of  the  diakopted 
antenna.  This  “complete  impedance  matrix”  relates  the  total 
scalar  potentials  to  all  the  impressed  currents. 

Interconnection  of  the  structure  elements,  which  requires 
equal  scalar  potentials  at  the  interconnected  terminals  and 
continuity  of  the  junction  currents,  is  formulated  by  an  inter¬ 
connection  matrix,  in  this  manner  a  system  of  linear  equa¬ 
tions  is  obtained  which  yields  the  junction  currents  and  the 
input  impedance  of  the  antenna. 

II.  CURRENT  COUPLING  BETWEEN 
STRUCTURE  ELEMENTS 

A.  Structure  Elements  with  One  Terminal 

Consider  one  of  the  capacitor  plates  of  the  dipole  in  Fig.  2 
separated  from  the  other  elements  and  suspended  in  space 


G(r.?')  = 


exp(-/*|f-r’l) 

1  —  » 


k  =  2irA 


where  f  is  the  position  vector  of  the  charges  and  currents 
on  the  surface  elements  dS,  and  f  is  the  point  of  observation. 

The  quantities  i(r)  and  q(f)  must  satisfy  the  following  two 
equations  on  the  surface  of  the  element  outside  the  contact 
area  o : 

E(?)  XdS  =  -[  juA(?)  +  V0(r)]  X  dS 

=  0  (boundary  condition)  (3) 

V  •  f(r)  +  /coq(f)  =  0  (continuity  condition).  (4)  j 


ine  condition  mat  current  fiux  through  the  boundary _ 

r  of  the  contact  area  o  is  the  continuation  of  the  impressed 
current  /  can  be  stated  as 


i 


i{r)-i(f)dr  =  l , 


(5) 


where  iff)  is  a  unit  vector  tangential  to  the  surface  S  and 
normal  to  r.  The  current  and  charge  distributions  1(f)  and 
q(f)  due  to  the  impressed  current  I  are  termed  as  “dominant” 
distributions  since  the  currents  due  to  field  coupling  between 
the  elements  are,  in  general,  relatively  small.  From  the  bound- 


GOUBAU  tt  alt  DIAKOPT1C  THEORY  FOR  MULTIELEMENT  ANTENNAS 


17 


ary  condition  (3) 


/  [/wI(f)+V«r))-J(?)dS  =  0. 
' s 


(6) 


The  surface  of  integration  S  is  the  surface  of  the  element  with 
the  exclusion  of  the  contact  area.  Using  the  relations 

V0  •  1(f)  =  V  •  (0(r)i(r))  -  0(r)V  •  I(r) 

=  V  •  (<Mf)l(}))  +  /«0(r)<?(r)  (7) 

and  applying  Gauss’s  theorem  one  obtains  from  (6) 


/«  /  t-4(r)  •  I(r)  +  0(r)q(')] 

■S 


-/ 


V  *  WfW))dS 


0(r)/(f)’  f(r)dr. 


(8) 


If  the  contact  area  o  is  sufficiently  small  0  can  be  considered 
constant  within  the  contact  area  and  on  T.  Thus  with  (5), 
(8)  reduces  to 

/«  f  l*r)  ’  i(r)  +  0(r)q(r")]  dS  =  <E»/  (9) 

Js 

where  0  is  the  scalar  potential  at  the  contact  area. 

The  ratio  between  4>  and  /  can  be  used  to  define  an  im¬ 
pedance  which  shall  be  termed  “intrinsic  impedance.”  If  A 
and  <t>  are  expressed  by  the  current  and  charge  distribution, 
the  intrinsic  impedance  of  the  element  is 

d> 

Z~- 

I 


/ 2 


*(O  +  0(')d(')]  dS 


4ir 


G(f,r') 


Hr)  *  Hr) 

l 2 


<f(r)q(r) 

k 2  Q 2 


dS’  dS 


(10) 


Both  terms  in  (1 1)  result  from  the  second  term  (charge  term) 
in  the  integrand  of  (10).  The  first  term  in  the  integrand 
(current  term)  yields  a  lowest  order  contribution  and, 
therefore,  vanishes  for  to  -*  0.  The  first  term  1/fcoC  in  (11)  is 
the  one  to  be  expected.  The  second  term  represents  a  negative 
resistance  of  -30  Q  for  a  conductor  in  free  space  (q  =  120sr) 
and  is  not  quite  obvious.  It  is  brought  about  by  the  fact  that 
an  impressed  current  produces  a  charge  on  the  element  with¬ 
out  a  countercharge,  in  contrast  to  a  Maxwell  source.  If  the 
scalar  potential  is  expanded  in  a  power  series  in  o>  one  obtains 


<Kr)~ 


-L / 

4ire 


G(r,  f')q(r)dS' 


1 

4  « 


dS'—jk 


I  q(r)dS'+ 
Js 


The  first  term  of  this  expansion  is  the  static  potential  of  the 
charges.  The  second  term  which  is  independent  of  r  repre¬ 
sents  a  potential,  termed  “background”  potential  0q,  which 
is  uniform  in  space  and  has  no  gradient.  This  means  that  0q 
does  not  produce  a  field.  It  is  this  background  potential  which 
produces  the  -30  12  term  in  (11).  When  the  element  which 
we  assumed  to  be  suspended  in  space  is  within  the  antenna 
structure,  the  background  potential  is  compensated  because 
the  total  of  the  combined  charges  on  all  the  other  elements  is 
equal,  but  opposite  in  sign,  to  the  charge  of  the  considered 
element.  The  background  potential  can  be  avoided  if  the  re¬ 
tarded  scalar  potential  is  redefined  as  a  modified  potential  0: 


*  1  I  •  , 

0  =  0-00= -  /  G(r,  r)q(r)dS' 

4  ire  Js 


(12) 


where 


G(r.  r)  =  -  -  -  +)k.  (13) 

I  r  —  r  | 

This  modified  scalar  potential  will  be  used  throughout  the 
paper.  It  is  preferred  over  the  conventional  potential  0  since  it 
will  lead  to  expressions  for  the  intrinsic  impedance  which  are 
closer  to  physical  expectation. 

Use  of  0  is  legitimate  as  it  does  not  conflict  with  Maxwell’s 
theory.  Since  V0  —  V0,  the  boundary  condition  (3)  and  the 
dominant  current  distribution  derived  from  it  remain  un¬ 
changed  if  0  is  replaced  by  0.  For  a  Maxwell  system,  0  and 
4  are  identical  since  fq  dS  extended  over  the  surface  of  the 
entire  structure  is  zero.  The  intrinsic  impedance  of  a  structure 
element  with  one  terminal  becomes 


where  Q  =  ///to  is  the  total  charge  on  the  element.  The  current 
and  charge  distribution  functions  I/I  and  q/Q  are  solely  deter¬ 
mined  by  the  geometry  of  the  element  and  the  location  of  the 
coupling  area. 

When  the  intrinsic  impedance  is  computed  with  (10)  for 
a  conductor  of  any  shape,  for  extremely  low  frequencies  it 
takes  the  form 


w-o  fuC 


(11) 


where  C  is  the  static  capacitance  of  the  element. 


M(f)  *  i(r)+  0(r)q<r)J  dS 

i(f) 


=  -  =  f 

/  i2  4 

1  A.-  d*r)q(?') "]  __ 

J  d*  dS- 


(14) 


Equation  (14)  represents  a  stationary  formulation  of  the 
intrinsic  impedance.  This  means  small  errors  in  the  dominant 
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current  distribution  have  only  a  second-order  effect  on  the 
intrinsic  impedance  (see  Appendix  I). 

Excitation  by  an  impressed  current  I  at  the  terminal  can  be 
considered  equivalent  to  excitation  by  an  oscillatint  charge 


/ 


which  is  placed  above  the  contact  area  at  a  distance  d  -*  0  as 
shown  in  Fig.  4.  The  charge  on  the  contact  area  o  consists 
essentially  of  the  image  charge  Q.  As  a  result  a  charge  +Q  is 
distributed  over  the  surface  area  of  the  structure  element  be¬ 
cause  the  net  charge  on  the  element  must  be  zero.  Under  cer¬ 
tain  conditions  charge  excitation  has  certain  practical  advan¬ 
tages  in  the  calculation  of  the  dominant  current  distribution. 

The  intrinsic  impedance  Z  of  an  element  with  one  terminal 
can  be  represented  by  a  lumped-element  circuit  as  shown  in 
Fig.  5.  For  low  frequencies,  i.e.,  when  the  dimensions  of  the 
element  are  small  compared  with  the  wavelength,  C  and  L  can 
be  considered  constant,  while  R  increases  proportionately  with 


Z  = 


- H  juL  +  R(co2). 

jioC 


(16) 


B.  Structure  Elements  with  Two  or  More  Terminals 

A  structure  element  with  two  terminals  such  as  the  cylin¬ 
drical  conductors  of  Fig.  2  has  two  dominant  current  distri¬ 
butions;  one  associated  with  each  of  the  independently  im¬ 
pressed  terminal  currents  (Fig.  6).  Each  dominant  current 
distribution  produces  a  scalar  potential  at  both  contact  areas. 
If  i  and  4>2  i  are ^  the  potentials  at  the  terminals  1  and  2 
due  to  /( ,  and  •i’u,  d>22  those  due  to/2,  then  the  relationship 
between  the  total  potentials  4>j  and  4>2  at  the  terminals  and 
the  impressed  currents  can  be  written  as 

4>,  =4>,,  +  4>i  2  =  Z(  ,/j  +  ZI2/2 

4*2  ~  <b2  i  +  ^*22  =  Z21/j  +  Z22/2.  (17) 

For  a  structure  element  with  M  terminals  the  relationship 
between  the  terminal  potentials  and  the  impressed  currents  is 
formulated  by  anftf  X  M  intrinsic-impedance  matrix; 

[*]  =  (*!(/]  (18) 


L 

AAAAA 


Fig.  S.  Low  frequency  equivalent  circuit  for  Nagle  terminal  atructuie 


cr, 


V' 


Fig.  6.  Structure  element  with  two  terminals. 


*11"  *12  *22'  *12 


Fig.  7.  Equivalent  circuit  lot  two-terminal  structure  element 

scalar  potentials,  respectively.  Equation  (20)  is  derived  in 
Appendix  II.  The  symmetry  of  the  intrinsic-impedance  matrix 
Zjk  =  Zk,  is  evident  from  (20).  In  Appendix  I  it  is  shown  that 
(20)  is  a  stationary  representation  of  the  impedance  matrix 
elements. 

A  lumped-element  equivalent  circuit  for  a  structure  dement 
with  two  terminals  is  shown  in  Fig.  7.  For  sufficiently  low 
frequencies  the  capacitors  and  inductors  can  be  considered 
constant;  while  the  resistors  increase  with  cj1  . 

A  power  series  expansion  of  the  right  side  of  (20)  in  terms 
of  /to  yields 

1  , 
zik  =  ——  +  /«£/*  +  */*(<*>  )  +  "• 


where 


I 


M*(')  *  ift)  +  <M^)q/(r)]  dS 


juji 

4ir 


'A 


dS'  dS , 


(19) 


(20) 


I/  =  lu>Qj  is  the  impressed  current  at  the  /th  terminal  (with 
lk  =  0)  and  lk  =  iuQk  is  the  impressed  current  at  the  fcth 
terminal  (with  //  =  0).  The  quantities  h,  qj,  A y,  and  lk, 
4 *,  Ak,  9k  *f*  the  corresponding  dominant  currents  distri¬ 
butions,  dominant  charge  distributions,  retarded  vector,  and 


where  C  is  independent  of  /,  k  (in  the  low  frequency  limit  the 
charge  distribution  on  a  given  structure  element  is  independent 
of  the  location  of  the  contact  area).  L)k  and  Rjk  arc  positive 
for  /  =  A,  but  in  general  are  negative  for  /  # ft.  In  this  case  the 
term  /col,*  may  be  counted  as  a  higher  order  contribution  to 
the  capacitive  term  1  //coC.  In  Fig.  7  this  assumption  has  been 
made  and  an  inductor  L\ 2 ,  which  should  appear  in  series  with 
the  capacitor  C,  has  been  omitted.  The  resistor  Rl2  of  this 
circuit  branch  may  be  negative,  but  its  magnitude  wil  be 
smaller  than  that  of  the  resistors  Rlt  -  Rn  and  R22  Rt2 
in  the  horizontal  branches  of  the  circuit. 

III.  FIELD  COUPLING  BETWEEN  STRUCTURAL 
ELEMENTS 

We  now  consider  a  diakopted  structure  and  arbitrary  cur¬ 
rents  impressed  at  the  terminals.  The  capadtivety  loaded 
dipole  of  Fig.  2  may  serve  as  an  example.  The  terminals  are 
identified  by  a  superscript  i  and  a  subscript  ft;  the  superscript 
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referring  to  the  number  of  the  element  and  the  subscript  re¬ 
ferring  to  the  number  of  the  terminal  on  the  element.  If  there 
were  no  field  coupling  between  the  elements,  the  current  dis¬ 
tributions  on  all  the  elements  would  be  the  dominant  distri¬ 
butions  associated  with  the  impressed  currents. 

The  field  of  a  dominant  current  distribution  is  non- 
Maxwellian  since  the  associated  net  charge  is  nonzero.  If  a 
current  /*'  is  impressed  at  the  terminal  (* ),  the  non-Maxwellian 
field  of  the  dominant  current  and  charge  distribution  i*  , 
qk  induces  currents  on  all  the  other  elements.  The  scatter 
fie'-H  excited  by  these  induced  currents  are  Maxwellian,  since 
induced  current  distributions  have  no  net  charge.  These 
"first-order''  scatter  fields  excite  second-order  scatter  fields 
and  so  on.  each  higher  order  having  a  greatly  reduced  ampli¬ 
tude.  All  these  scatter  fields  when  summed  up  form  a  multiple- 
scatter  field  which  is  Maxwellian.  The  currents  and  charges 
associated  with  the  multipie-scatter  field  are  distributed  over 
all  the  surfaces  S"  (including  S1)  and  shall  be  denoted  Si*1", 
Sqkln ,  the  superscripts  and  subscripts  indicate  that  they  are 
produced  by  the  impressed  current  /*'  and  located  on  the 
element  n.  In  calculating  6/*'"  and  5<7*‘"  the  structure  ele¬ 
ments  are  assumed  open-circuited  at  their  terminals.  Hence 
the  total  flux  of  the  scatter  currents  through  the  rim  of  the 
contact  areas  is  zero.  Scatter  currents  do  not  contribute  to  the 
junction  currents. 

The  total  field  generated  by  /*'  (all  other  terminals  open- 
circuited)  satisfies  on  every  element  the  boundary  conditions: 

(M/l*'  +  bAk‘)  +  7(«*'  +  60*'))  X  dSn  =  0, 

n  =  1,  i,  •••,  N  (21) 

where  Ak‘,  0*'  are  the  retarded  potentials  of  the  dominant 
current  and  charge  distribution  i*',  <?*';  and  bAk‘,  80*'  are 
thoee  of  the  scatter  current  and  charge  distributions  6 »*'”, 
bqk"  combined.  N  is  the  number  of  elements. 

Since  the  electric  field  of  i*',  qk  satisfies  the  boundary 
condition  on  s',  it  follows  from  (21)  for  n  =  i  that 

(JubAk‘  +  V60*')Xd?  =  O. 


Furthermore,  from  the  boundary  conditions  (21),  using  the 
relations  (7)  and  Gauss’s  theorem,  follows 

/  [(5*'  +  6Ak‘)  •  blk,n  +  (0*'  +  60*')5<j*'"J  dSn  -  0, 

Jsn 

for  every  n  including  i.  (25) 

The  right  side  of  (25)  is  zero  since  for  scatter  currents  the 
rim  integral  of  Gauss’s  theorem  in  (5)  is  zero  (see  Appendix  IV). 

With  (23),  (24),  and  (25)  one  obtains  the  “backscatter” 
impedances 

7  t,i  4>*  *'’'(F) 

*k,k  (F)=  — — — — 

h 

=  -/«(- l  /  («X*'  •  bikin 

V*  /  n=  1  JSn 

+  60* '6  <?*'")  dS",  (26) 

which  has  to  be  added  to  the  diagonal  terms  of  the  intrinsic- 
impedance  matrix  Zk,k%i ,  using  the  notation  of  this  section. 
Generalization  of  (26)  to  obtain  the  scatter  field  contributions 
to  the  off-diagonal  terms  is  straightforward  and  yields 


Z*/'(F)  = 


*kj\F) 


n  =  1,  -,i,  -,N  (21)  '"/w-ttS  /  0Ak‘-&I/n 

—  ,  *  ,  / *  /;  n=  i  gw 

where  Ak,  0*'  are  the  retarded  potentials  of  the  dominant 

current  and  charge  distribution  /*',  <?*';  and  6 Ak  ,  80*'  are  +60* '6 q/n)dS”.  (27) 

thoee  of  the  scatter  current  and  charge  distributions  6/*'”, 

bqk"  combined.  N  is  the  number  of  elements.  For  k  =/  (27)  transforms  into  (26). 

Since  the  electric  field  of  i*',  qk  satisfies  the  boundary  Let  us  now  determine  the  potential  4>*>m  ,l(F)  produced 
condition  on  S1,  it  follows  from  (2 1 )  for  n  =  i  that  at  the  terminal  (*)  by  the  current  Ime  impressed  at  the  terminal 

.  (  ').  From  the  boundary  condition  (21) 

(ju>bAk  +  760*')  X  dS1  —  0. 

Thu*'  /  [MAm-'  +  bAj)  +  V(0m'  +  60m')]  •  VdS*  =  0, 

JS‘ 

I  [(/cj6A*'  +  780*')  •  f*'|  dS*  —  0,  (28) 

Js‘ 

i  =  1 ,  — ,  Af;  k  =  1 ,  — ,  Af,.  (22)  and 

Using  the  relations  (7)  and  Gauss’s  theorem  one  obtains  the  r  _  _  „ 

"backscatter”  potential  due  to  the  field  interaction  of  the  j  [/w(/f*' +  8/1*')  +  7(0*' +  60*'))  *  bimln  dS"  =  0, 
excited  element  with  the  other  elements:  ■'£" 

♦»/'(^)/*'*/«  /  («/<*'•  y  +  «0*V)<«'.  (23)  U<’* 

s*  where  the  second  equation  holds  for  every  n  including  i.  The 

The  letter  F  indicates  field  coupling;  the  first  pair  of  indices  potentials  bAk  and  60*'  characterize  the  scatter  field  which 
(*)  refers  to  the  terminal  at  which  <!>  is  determined,  and  the  would  be  excited  by/*'. 

second  pair  to  the  terminal  of  the  impressed  current  which  As  before,  we  apply  relations  (7)  and  Gauss's  theorem,  and 
produces  this  potential.  obtain  from  (28) 

As  shown  in  Appendix  111 


/  («/(*'•  I*' +  60*'**')  rfS' 

V 

-2  /  <x*' •  *y" + 0* wvs" 

»>l 


KmU(F)lk‘  =  /«[  fsym'  •  y  +  ♦e.'ffcW 

+  |(6Am'*y  +  6^V)^'j  (30) 
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and  from  (29) 

0=  (  {Ak‘  •  Simln  +  0*,5<jm,")</Sn 

+  /  ( &Ak‘-6lm"'  +  6ik,6qm‘n)dSn .  (31) 

The  first  term  in  (30)  represents  the  contribution  to  the  ter¬ 
minal  potential  from  the  non-Maxwellian  field  of 

the  dominant  current  and  charge  distribution  iml,  qm‘,  and  the 
second  term  that  from  the  scatter  current  and  charge  distri¬ 
butions  b~imln  ,bqmln . 

As  shown  in  Appendix  III 

/(6Am'-iV  +  «0mV)rfS' 

JS‘ 

=  2  /  (Ak,-S!m"'+k‘Sqm,n)dSn.  (32) 

Expressing  4>*,m,’,(F)  in  terms  of  an  impedance 

*k,mi'\F)  =  Zk,m\F)lJ,  (33) 

and  using  (30),  (31),  and  (32),  the  field-coupling  impedance 
between  the  terminals  (*)  and  (J)  becomes 

Zk,m,J(F)  =  777T  /  &•*  ’  *k‘ +  im'qAdS1 

‘k  'm  L  JS 1 

~  2  /  <***'  •  6im'n+bk,bqm'n)dSm  , 

n=  1  S"  J 

i*l.  (34) 

Equations  (27)  and  (34)  yield  the  elements  of  the  field-cou¬ 

pling  impedance  matrix  [2(F)] ,  which  relates  the  scalar  poten¬ 
tials  4>k‘(F)  at  the  terminals,  caused  by  field  coupling,  to  the 
impressed  currents: 

wni  =iz(F)M, 

Win* 2  2  Km  \n 

1=1  «n=l 
N  Ml 

=  22  (35) 

/■  l  m«i 

IV.  COMPLETE  IMPEDANCE  MATRIX  OF  THE 
DIAKOPTED  ANTENNA 

The  intrinsic-impedance  matrices  of  the  individual  struc¬ 
ture  elements  can  be  combined  into  a  diagonal  block  imped¬ 
ance  matrix  (2(C))  by  writing  the  matrix  elements  Zkj  (20) 
in  the  form  Zkj‘(Q.  The  superscripts  /  identify  the  terminals 
k  and  /  as  belonging  to  the  element  1;  the  letter  C  indicates 
current  coupling.  The  block  matrix  (2(C)),  whose  elements 
2*  jU(C)  are  zero  for  /  #  /,  is  the  “current-coupling  matrix’’ 


of  the  diakopted  gystem  and  relates  the  terminal  potentials 

M,  * 

**'(o=2  2*/'(oV 


due  to  current  coupling  to  the  M /  impressed  currents  of  the 
element  i. 

The  sum  of  the  matrices  [2(C)]  and  [2(F)] ,  i.e., 

[2]  =[2(C)] +[Z(/0]  (36) 

forms  the  “complete  impedance  matrix’’  of  the  diakopted 
antenna,  which  represents  the  relationship  between  the  total 
terminal  potentials 


=  V(0  +  4>*'(F)  =  2  2  Km1,1 

1=1  m  =  1 

produced  by  current  and  field  coupling  to  all  impressed  cur¬ 
rents.  In  matrix  form 


[4»]  =  [^][/]. 


(37) 


If  the  matrix  elements  ZkJ*{C)  (20)  and  Zkj\F)  (27)  are 
added,  the  resulting  elements  2(kl/<■^  have  the  same  form  as 
those  which  pertain  to  field  coupling  between  different  struc¬ 
ture  elements  (34).  In  other  words  (34)  can  be  used  as  general 
formulation  for  all  elements  of  the  complete  impedance  ma¬ 
trix  of  the  diakopted  system,  and  the  condition  i  *  I  can  be 
dropped. 

The  symmetry  of  the  [2]  matrix,  i.e., 


7  U  —  7 
Lk,m  m.k 


(38) 


can  be  easily  verified,  by  expressing  in  (34)  the  vector  and 
scalar  potentials  by  the  current  and  charge  distributions  ac¬ 
cording  to  ( 1 )  and  ( 1 2). 

Equation  (34)  represents  a  stationary  formulation  of  the 
matrix  elements  of  [2],  This  means  first-order  errors  in  the 
current  and  charge  distributions  lead  to  only  second-order 
errors  in  the  impedances  (Appendix  I). 

Calculation  of  the  radiation-coupling  impedances  according 
to  (27)  and  (34)  requires,  in  principle,  computation  of  the 
scatter  current  and  charge  distributions.  However  numerical 
results  obtained  with  this  theory  indicate  that  coupling  by  the 
scatter  currents  is  a  negligible  effect.  It  has  been  found  that 
coupling  by  the  junction  currents  prevails  over  field  coupling, 
and  field  coupling  by  the  non-Maxwellian  fields  (radiated  by 
the  dominant  current  distribution)  dominates  over  that  by 
the  Maxwellian  fields  (generated  by  the  scatter  currents).  In 
principle  field  coupling  effects  by  the  scatter  currents  can  be 
obtained  with  an  iterative  procedure,  which  is  not  discussed 
here. 

If  coupling  by  the  Maxwellian  scatter  fields  is  neglected, 
the  formula  for  the  elements  of  the  complete-impedance 
matrix  for  the  diakopted  system  reduces  to 

Z*.»'''“77rT  /  Mm'-ia'  +  fm'daW  (39) 
VV  JS< 

This  means  ail  the  matrix  elements  can  be  computed  from  the 
dominant  current  distribution.  It  is  worth  mentioning  that 
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formulation  of  the  antenna  problem  in  terms  of  this  approx¬ 
imate  impedance  matrix  can  also  be  derives  in  a  more  con¬ 
ventional  manner  by  an  application  of  the  method  of  momenta 
(3].  (41.  For  this  purpose  one  would  formulate  the  electric 
field  integral  equation  of  the  entire  (interconnected)  antenna 
structure  and  reduce  this  equation  to  a  linear  system  by  em¬ 
ploying  the  Galerkin  version  of  the  method  of  moments, 
while  choosing  the  dominant  current  distributions  of  the  vari¬ 
ous  structure  elements  as  (subsectional)  basis  and  weight  func¬ 
tions.  After  transformation  of  the  linear  system,  by  applica¬ 
tion  of  (7),  into  a  matrix  equation  between  the  potentials  and 
currents  at  the  contact  areas,  the  matrix  coefficients  are  found 
to  be  identical  with  (39).  The  choice  of  the  dominant  current 
distributions  as  basis  functions  will  alleviate  the  necessity  for 
subdividing  a  given  antenna  into  many  small  segments.  Since 
each  of  these  current  distributions  satisfies  the  conditioning 
ftan  =  0  at  its  self-element  (and  since  current  coupling  domi¬ 
nates)  it  can  be  expected  that  structure  elements  of  compara¬ 
tively  large  size  will  already  yield  accurate  results.  In  other 
words  the  accent  of  the  problem  is  shifted  from  salving  an 
integral  equation  for  the  entire  antenna  to  determination  of 
the  dominant  current  distributions  of  its  structure  elements 
which,  in  many  cases,  can  be  reduced  to  canonical  problems. 
Moreover,  because  of  the  stationary  properties  of  the  matrix 
coefficients  (39),  the  dominant  current  distributions  have  to 
be  known  only  approximately  in  order  to  obtain  impedance 
values  of  good  accuracy.  The  numerical  example  discussed  in 
Section  VI  confirms  these  predictions. 

V.  INTERCONNECTION  OF  THE  STRUCTURE  ELEMENTS 
AND  INTERCONNECTION  MATRIX 


-4  «Lg  l *  !*  »’  > 

i4.  i?  U  H  tt  ii 

Fig.  8.  Thin  wire  dipole  treated  as  diakopted  four-clement  system. 


relationships  are 

(*]  =  (*](/] 

(diakopted  antenna] 

(40) 

[*r=(z]vr 

(assembled  actual  antenna] . 

(41) 

Primed  quantities  refer  to  currents  and  potentials  at  various 
interconnections  of  the  actual  antenna.  Requirements  a),  b), 
and  c)  represent  Kirch  offs  laws  for  interconnected  structures 


and  can  be  written  as 

C  i]-cci(/]' 

(42) 

(*]'  =  (  <:],[<*>] 

(43) 

(<*>]/[ /]’  =  [*]«{ /]. 

(44) 

[C)f  represents  the  transpose  of  [C],  Note  that  <£'  denotes 
potential  differences,  i.e.,  voltages. 

From  (42),  (43),  and  (44)  the  impedance  of  the  actual 
interconnected  structure  can  be  written  as 


The  requirement  for  the  diakopted  structure  with  im-  —  [C] .  (45) 

pressed  currents  to  be  identical  in  performance  with  the 

assembled  antenna  are  as  follows.  The  following  example  shows  how  [C]  and  [Z] '  are  obtained. 


a)  The  sum  of  the  impressed  currents  is  zero  at  every  junc¬ 
tion  between  the  structure  elements.  This  requirement 
assures  that  the  Kirchoff  condition  of  current  continuity 
is  satisfied  and  that  the  field  of  the  assembled  antenna 
is  Maxwellian. 

b)  The  scalar  potentials  at  interconnected  terminals  are 
equal. 

c)  The  potential  difference  between  the  input  terminals  is 
equated  with  the  driving  voltage  of  the  antenna. 

It  should  be  emphasized  that  the  continuity  of  current  at  a 
junction  needs  to  be  imposed  only  on  the  impressed  cur¬ 
rents  which  are  in  turn  the  fluxes  of  the  dominant  current 
distributions  through  the  contact  areas.  The  scatter  currents 
do  not  contribute  to  the  terminal  currents  since  they  are 
Maxwellian  and  have  no  net  charge.  Thus  the  satisfaction  of 
boundary  conditions  (21)  is  consistent  with  the  satisfaction  of 
terminal  condition. 

Imposing  these  junction  conditions  the  matrix  equation 
(36)  yields  a  system  of  linear  equations  for  the  unknown  junc¬ 
tion  currents  and  the  input  impedance  of  the  antenna.  Using 
network  theory  concepts  the  reduction  of  (37)  to  this  linear 
system  of  equations  by  enforcing  the  junction  conditions  can 
be  formulated  with  a  connection  matrix  [C]  which  reduces 
the  number  of  potentials  and  currents  of  the  diakopted  struc¬ 
ture  to  those  of  the  actual  structure  [5] .  The  potential-current 


VI.  EXAMPLE 

As  an  example  we  apply  the  diakoptic  theory  to  an  ordi¬ 
nary  thin-wire  dipole  antenna  and  compare  the  results  with  the 
data  available  in  the  literature.  To  obtain  a  multielement 
structure  we  cut  each  wire  in  half,  as  shown  in  Fig.  8,  and  con¬ 
sider  each  half  as  a  structure  element.  The  diakopted  dipole  is 
thus  modeled  by  two  structure  elements  with  one  terminal  and 
two  structure  elements  with  two  terminals,  so  that  the  total 
number  of  terminals  is  six.  The  complete  impedance  matrix 
of  the  diakopted  structure  [Z]  is  therefore  a  6  X  6  matrix. 
However  there  are  only  eight  different  impedances  because  the 
four  structure  elements  have  been  assumed  to  be  alike. 

Using  the  enumerations  of  Fig.  8  the  matrix  equation  (37) 
has  the  form 


*.3 


4>j* 


V 
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with 


7  —  7  33  _  7  11  7  11  —  7  22  —  7  22  7  44 

*0”*il  “*22  ”*11  *22  *11  ”*22 

7  —  7  31  —7  13  —  7  12  —  7  21  —  7  24  7  42 

*1”*12  ”*21  ”*12  “*21  ”*I2  ”*21 

7  —  7  11  —  7  11—7  22  —  7  22 

*2”*21  ”*12  ”*21  *12 

7  —  7  31  —  7  13  —  7  12  —  7  21  —  7  12  —  7  21—7’  24  —  7  42 

*3”*11  ”*11  ”*11  *11  ”*22  ”*22  ”*22  ~  *22 

y  —  7  32—7  23  —  7  14  —  7  41 

*4  ”*12  ”*21  “*12  *21 

7  —7  12  —  7  21 

*5  ”  *21  ”  *12 

7  —-7  32  —  7  23  —  7  14  —  7  41 

*6  *11  ”*11  ”*22  *22 

7  —7  34  —  7  43 

*7  ”  *12  ”  *21 


The  darkened  portion  of  the  impedance  matrix  Z  is  the 
current-coupling  matrix  Z(C). 

The  interconnection  conditions  require 


ample  are  the  same  for  all  the  elements,  can  be  approx¬ 
imated  by  linear  current  distributions  (uniform  charge 

distribution). 


/, 3  =-/,«=/,  $,3=4)21=4,, 

I72  =-/, 1  =  -/0  -4>22  =  K0 

/24  =-/,*=/ J  4»24  =  4>,2  =*, 

where  /0,  K0  are  input  current  and  driving  voltage  of  the 
antenna. 

Because  of  the  symmetry  of  the  antenna 
/,  =  -~/2 ;  =-4>2;4>2J  = 

Current  and  voltage  matrix  of  the  interconnected  antenna  are 

1/1  = 

Thus  the  interconnection  matrix  becomes 


and  the  impedance  matrix  of  the  assembled  antenna  is 


Although  the  latter  approximation  is  rather  crude,  one 
should  expect  reasonable  results  if  the  wire  sections  are  short 
compared  with  the  wavelength  because  all  the  impedance 
formulas  are  stationary  expressions.  Linear  current  distribu¬ 
tion  permits  analytic  formulations  of  all  the  impedances 
Z0,  Z ! ,  Z2,  etc.,  and  numerical  calculations  with  a  pocket 
calculator  (such  as  HP  25).  The  results  obtained  are  presented 
in  Fig.  9.  The  curves  are  plots  (from  a  table  by  King  [6])  of 
the  real  and  the  imaginary  part  of  the  input  impedance  of  a 
dipole  with  In  (2//p)  =  5  as  a  function  of  (3  =  kl;  21  is  the  total 
length  of  the  dipole  and  p  is  the  wire  radius.  The  crosses  mark 
the  values  of  the  input  impedance  as  obtained  from  (46)  with 
the  above  assumptions.  For  0  =  kl  <  0.8  the  deviation  of  the 
real  part  of  the  input  impedance  from  King’s  data  is  less  than 
10  percent,  and  for  the  imaginary  part  deviation  is  less  than 
1  percent.  From  this  one  can  conclude  that  the  linear  approx¬ 
imation  for  the  dominant  current  distribution  is  adequate  if 
the  length  of  a  wire  section  is  <1/15  A.  This  has  been  con¬ 
firmed  by  computer  results  obtained  with  each  dipole  leg 
diakopted  into  four  equal  sections.  These  results  are  marked 
in  Fig.  9  by  dots  and  are  in  good  agreement  with  King’s 
(curves)  even  beyond  the  second  resonance  of  the  antenna. 

Fig.  10  gives  an  indication  of  the  convergence  of  numerical 
results  as  segmentation  of  a  given  antenna  increases.  A  dipole 
antenna  is  assumed  with  each  leg  diakopted  into  N  equal  parts; 
the  relative  impedance  errors  \AR/R  I  and  \AX/X  I  are  plotted 
versus  the  antenna  parameters  SI  =  2  In  (2 I/p)  for  0  =  2  and 


Hence 


(Z0-Z1) 

(2Zj-Z2-Z4) 

(2/3  —  z2  —  Z4) 

(2Z0-2Zt  ~Zs  +  2Z6  —  Z7) 

- Vhzlizvt 

l0  2Z0  — 2Z,  -ZS  +  2Z6-Z7 


(46) 


For  the  numerical  calculation  of  the  impedances  the  following 
amplifying  assumptions  have  been  made. 

a)  Coupling  by  scatter  currents  is  negligible. 

b)  The  dominant  current  diatributiona,  which  in  this  ex¬ 


various  N.  The  impedance  errors  were  determined  by  compar¬ 
ison  with  King’s  results  [6] .  In  general  these  errors  decrease 
with  increasing  N.  But  the  behavior  of  the  curves  in  the  region 
of  smaller  fl(=10  •••  12)  suggests  that  for  any  12  there  is  an 
optimum  value  of  /V  for  which  the  errors  are  at  a  minimum. 


V 
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Fig.  9.  Comparison  of  dipole  impedance  calculated  with  diikoptic 
theory  versus  King. 


Fig.  10.  Difference  between  King's  results  and  dialcoptic  theory’s 
results. 

while  •  further  increase  in  N  would  tend  to  reduce  accuracy. 
As  one  should  expect  the  optimum  N  increases  with  f2,  i.e., 
as  the  antenna  becomes  more  slender.  The  length  of  the  wire 
sections  can  be  substantially  increased  if  their  dominant  cur* 
rent  distributions  are  approximated  by  sinusoidal  rather  than 
linear  functions.  Fig.  11  shows  the  input  admittance  of  a 
folded  dipole.  The  solid  curves  were  taken  from  a  table  by 
King  and  Harrison  (7).  The  present  method  was  evaluated 
counting  the  horizontal  members  of  the  antenna  as  one  struc* 


Fig.  11.  Comparison  of  folded  dipole  admittance  calculated  w:th  dia¬ 
koptic  theory  versus  King  and  Harrison. 


ture  element  each,  while  each  vertical  member  was  diakopted 
into  two  sections.  For  the  dominant  current  distributions  on 
these  sections  the  sinusoidal  approximation  was  used.  Despite 
the  comparatively  large  size  of  the  vertical  structure  elements, 
computed  impedance  values  were  found  to  be  in  good  accu¬ 
racy  up  to  and  beyond  the  second  resonance. 

VII.  CONCLUSION 

The  major  advantage  of  the  diakoptic  theory  for  multi¬ 
element  antennas,  as  discussed  in  this  paper,  is  that  the  prob¬ 
lem  of  determining  the  current  distribution  on  the  antenna 
need  not  be  solved  for  the  structure  as  a  whole  but  only  for 
the  individual  structure  elements.  Excitation  of  each  structure 
element  is  ascribed  to  the  currents  at  its  junction  with  adjacent 
elements  and  to  the  fields  of  the  surface  currents  on  all  the 
other  elements.  The  current  distributions  produced  by  the 
junction  currents  have  been  termed  dominant  current  distri¬ 
butions  because  they  constitute  the  major  portion  of  the 
currents  on  the  composite  antenna  structure.  The  remainder 
of  the  currents  are  made  up  by  scatter  currents  which  are 
produced  by  field  coupling.  Field  coupling,  in  first  approx¬ 
imation,  is  determined  by  the  dominant  current  distributions; 
while  coupling  by  the  scatter  currents  in  general  is  negligible. 
Introduction  of  impedances  for  the  characterization  of  struc¬ 
ture  elements  and  their  interaction  permits  use  of  network 
theory  concepts  for  the  determination  of  the  junction  currents 
and  the  input  impedance  of  the  antenna.  Formulation  of  all 
impedances  by  stationary  expressions  renders  the  results 
insensitive  to  computational  errors  in  the  current  distributions. 
As  demonstrated  by  the  example  given  in  the  paper,  even 
rather  crude  approximations  to  the  dominant  current  distri¬ 
butions  can  yield  goocj  results. 

APPENDIX  I 

Proof  for  the  Stationary  Formulation  of  the  Impedances 

1)  To  prove  that 
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(A  •  i  +  jq)  dS 


(Al) 
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represents  a  stationary  formulation  of  the  intrinsic  impedance,  or 

we  assume  that  the  dominant  current  distribution  1  has  an 

error  of  Ai.  The  corresponding  errors  of  <j,  X,  and  4  (hall  be  f  _ 

denoted  Aq,  AA,  and  A$.  Then  /<•>  ]  (Ak  •  Alj  +  pkqj)  dS  -  0. 

S 

Z  +  AZ=^  I  1(1+  AA)  •  (f  +  Ai)  Furthermore 

1  5  / 

+  (*  +  A*X«  +  (A2)  J(AAk-ij  +  A4kqt)  dS 


+  (d  +  A*X<?  +  ^?)1  dS.  (A2) 

The  boundary  condition  for  the  correct  dominant  current 
distribution  yields 

f  [<joA  +  V0)  •  Ai]  <fS  =  0.  (A3) 


=  /  (A  i  •  Ai] 


k  +  4>i&qk)dS  -  o. 


From  (A10)  and  (All),  (A9)  reduces  to 


Since  the  dominant  current  distribution  is  the  continuation  ,  r 

of  the  impressed  current  which  is  assumed  to  be  unchanged,  A Z,k  = -  /  (A Ak  •  Alt  +  A4>kAqj)dS. 

A  i  (I.  at  ika  ta«MM4Ma1  am<1  /  A  3  \  ahm  Wa  ...aIMaa  Iaa  SLa  f  j  /y 


Ai  is  zero  at  the  terminal,  and  (A3)  can  be  written  in  the  form 


iu  I  (A  *  Ai  +  p 
•v 


AAq)dS  =  0. 


Using  the  relations 


I  A  •  AldS-  I  AA  •  idS 
•s 

/  0Aj  dS  =  /  A£/  dS 
■'s 


along  with  (A4)  one  obtains 


Thus  A  Zjk  is  of  second  order. 

(A4)  3)  To  prove  that  (34)  is  a  stationary  expression  for  the 

field -coupling  impedances  we  treat  the  assembly  of  discon¬ 
nected  structure  elements  like  a  single  body.  This  means  when 
a  current  is  impressed  on  terminal  (k)  we  consider  the  dom¬ 
inant  current  distribution  ik  together  with  the  associated 
scatter  currents  6ik  which  are  distributed  over  all  the  elements 
as  a  dominant  current  distribution  of  the  system.  The  coupling 
impedances  between  any  two  terminals  can  then  be  formu- 

(AS)  lated  like  mutual-intrinsic  impedances  (A8), 


I  (AA  •  f  +  A4>q)  dS  —  £(A'Ai+*Aq)dS  =  0.  (A6) 

Thus  from  (AI),  (A3),  and  (A6) 

AZ=—.  I  (AA  •  A  7  +  A$Aq)dS.  (A7) 

t  JS 

This  means  A  Z  is  of  the  second  order. 

2)  In  the  case  of  a  mutual-intrinsic  impedance 

ju  f  -  .  . 

Z*~77  /  (A*  *  if  +  9*4/)  dS  (A8) 

•kb  5 

both  the  dominant  current  distributions  ik  and  i,  may  have 
errors  Ai*  and  Ai/.  Thus 


zk.mi'\F)-~T7~T  I  l(Am'  +  6Am') 
him 


/W  f 

***=77  1 

hb  Js 


him  JZSn 

•(i*,  +  fifV)+(0m,  +  6«mf) 

*  (q k  +«<?*')]  dS.  (A13) 

The  error  A Zkma  produced  by  errors  in  the  current  distri¬ 
butions  ik‘,  6ik  and  im\ Siml  is  obtained  from  (A12), 

AZ.J-'CF)^-—  f  [(AAk‘  +  A£Ak‘) 
him  JXS" 

*  +  Afi'm1)  +  (Ad**  +  ASd*1) 

(A14) 

and  is  of  second  order.  This  relation  can  also  be  derived  from 
(34),  but  only  in  a  rather  cumbersome  manner. 


[(A4k  •  1/  +  A 4kqj) 
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+  (Ak  •  Ai,  +  4kAq,)  +  (Xk  •  Ai,  +A <fikAq,) J  dS.  Derivation  of  (19) 


Consider  a  structure  element  with  several  terminals  and 
let  *  and  /  be  any  two  terminals  where  currents  /*  and  1,  are 
Because  the  correct  dominant  .current  distributions  satisfy  the  impressed.  The  corresponding  dominant  current  and  chaise 
boundary  condition  EXdS  =  Q,  distributions  f*,  qk  and  i,,  q,  produce  the  fields  Ek  and  E, 

/which  utisfy  the  boundary  conditions 
(/«J*  +  $d*)  V;<fS-  0 

-  £*Xd5«-(/wAk  +  7#k)XrfS«  0  (A1S) 
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(A1 6)  Multiplying  (A23)  with  ik(r)  and  integrating  over  S* 


EjXdS  =  ~  (JuiAj  +  Vfy)  X  dS  =  0.  (A16)  Multiplying  (A23)  with  i 

Since  the  currents  ik  and  1/  are  tangential  to  the  surfaces  it  j  .  j^^dS* 

fallows  from  the  boundary  conditions  'si 

f  (j<Mk  +  Wk)-IkdS  =  0  (A  17)  =~ 

f  =~y  (siini 

I  (/wAt  +  V0fc)  •  If  dS  =  0.  (A  18)  4ir  Jsn  k 

'c 


Using  the  vector  identity 


V  *  <£)  =  7*  •  i +  0(V  •  i\ 

with  V  •  /  =  —/ox? 


=  Ji  fik‘(r)-Y  {  f>htn{?yG{i,r)dSn<f)dSi(r) 
4tr  it  jjr,  Js* 

=  —  2  /  6fW)  /  ?*)  dS‘(r)  dS"(r) 
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N  r 

=  2  /  5,V"(')  *  Ak(r)df‘; 


««I  , sn 

G(n?')=C(r,r). 

(A1 9)  Similarily  it  can  be  shown  that 


and  applying  Gauss's  theorem  as  in  (8)  of  Section  11,  (A17) 
and  (A18)  can  be  written  in  the  form 

—  I  [  V  *  (4>k‘k)]  dS  =  4>kklk 

JS 

=  ;'w  f  (Ak  *  ifc  +  $kqk)dS  (A20) 


-  A  V  •  (0Jc'/)]  -  4/fc/y 

JS 

=  /W  /  (Ak'Ij  +  4>kqi)dS,  (A21) 


where  4>ka  «  the  potential  at  the  terminal  k  due  to  lk  and  <?* 
and  that  at  the  terminal  j  due  to  ik  and  qk  Of  course  i* 
and  qk  are  induced  on  the  structure  by  Ik  at  the  fcth  terminal 
with  all  the  other  terminals  open-circuited,  i.e.,  with  other 
impressed  currents  set  equal  to  zero.  For  /  =  k,  (A21)  trans¬ 
forms  into  (A20). 

With 

*ik  =  Zjklk  (A22) 

one  obtains  from  (A21)  the  expression  for  Z;*  given  in  (19). 
The  expression  in  (20)  is  obtained  if  the  potentials  A  k  and  4>k 
are  expressed  by  (1)  and  (12),  respectively. 


APPENDIX  III 

Derivation  of  (23) 

Let  ik  be  a  dominant  current  distribution  on  the  surface 
S1  and  Sr1*'"  be  the  scatter  current  distribution  on  5"  produced 

by  ?*'(»=  1,2,  — ,  N). 

u  N  f 

6V(')  =  ~  2  /  5ikln(f)G(f,-r')dSn(f).  (A23) 
4ff  * jn 


/  8^VdS'=2  /  6 ntnk‘dSn. 

_ .  'cB 


sl  rt=  1  sn 

The  proof  for  (32)  in  the  body  of  the  paper  follows  the 
same  outline  as  given  above. 

APPENDIX  IV 

If  a  current  is  impressed  on  any  terminal  of  a  diakopted 
structure  there  will  be  capacitive  currents  between  the  contact 
areas  of  the  disconnected  elements  which  have  not  been  con¬ 
sidered  in  the  derivation  of  the  field  coupling  impedances. 
One  might  therefore  conclude  that  the  formulas  are  approx¬ 
imations  which  require  the  gaps  between  adjacent  contact 
areas  to  be  so  large  that  capacitive  currents  are  negligible. 
The  purpose  of  this  appendix  is  to  show  that  the  expressions 
for  Zk%ki,l{F)  and  Zktm'J(F)  are  correct  even  if  the  gaps  are 
infinitely  small. 

Fig.  12  shows  two  structure  elements:  a  cylindrical  rod  1 
and  a  disc  2  with  the  opposing  contact  areas  O21  and  a(  2 .  If 
a  current  is  impressed  on  the  terminal  (} )  of  the  rod  there  will 
be  a  potential  difference  between  o2 1  and  Oi 2 ;  which,  in  turn, 
produces  a  displacement  current  between  these  terminals.  The 
potential  difference  which  is  the  line  integral  of  the  electric 
potential  field  between  o2 1  and  ot  2  is  essentially  determined 
by  the  charges  on  the  contact  areas.  If  the  gap  is  made  smaller 
and  smaller  the  potential  difference  approaches  zero  and  the 
total  current  distribution  becomes  the  dominant  current 
distribution  of  the  interconnected  elements.  Because  of  the 
equivalence  of  current  and  charge  excitation  displacement 
currents  at  contact  areas  are  equivalent  to  impressed  currents. 
Thus  the  situation  discussed  above  is  the  excitation  of  a 
diakopted  structure,  not  by  one,  but  by  three  impressed  cur¬ 
rents.  To  produce  excitation  by  one  impressed  current  in 
accordance  with  our  theory,  the  displacement  currents  must 
be  compensated  so  that  there  is  no  current  flux  from  the 
contact  area  onto  the  surface  S  of  the  element  (S  by  definition 
does  not  contain  the  contact  areas  of  the  element).  The  mag¬ 
nitude  of  these  compensating  currents  does  not  enter  into  the 
analysis  because  if  the  impressed  currents  of  the  diakopted 
structure  are  identical  with  the  junction  currents  of  the  inter- 
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Fig.  1 2.  Compensation  of  capacitive  currents  at  contact  areas. 

connected  structure  there  are  no  displacement  currents  be¬ 
tween  adjacent  contact  areas  and  the  sum  of  all  the  com¬ 
pensating  currents  is  zero. 
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